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Abstract 

Recent results demonstrating the chaotic behavior of geodesies in non-homogeneous 
vacuum pp-wave solutions are generalized. Here we concentrate on motion in non- 
homogeneous sandwich pp -waves and show that chaos smears as the duration of these 
gravitational waves is reduced. As the number of radial bounces of any geodesic de- 
creases, the outcome channels to infinity become fuzzy, and thus the fractal structure of 
the initial conditions characterizing chaos is cut at lower and lower levels. In the limit 
of impulsive waves, the motion is fully non-chaotic. This is proved by presenting the 
geodesies in a simple explicit form which permits a physical interpretation, and demon- 
strates the focusing effect. It is shown that a circle of test particles is deformed by the 
impulse into a family of closed hypotrochoidal curves in the transversal plane. These are 
deformed in the longitudinal direction in such a way that a specific closed caustic surface 
is formed. 
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1 Introduction 



The widely known class of plane- fronted gravitational waves with parallel rays (pp -waves) has 
become a paradigm of exact radiative spacetimes in general relativity. Found by Brinkmann 
in 1923 [|l| and later rediscovered independently by several authors, it has been studied for 
decades. The metric of vacuum pp -waves can be written in the standard form 

ds 2 = 2 dCdC - 2 dudv - (/ + /) du 2 , (1) 

where f(u, () is an arbitrary function of the retarded time u and the complex coordinate 
( spanning the plane wave surfaces u = Uq = const. When / is linear in (, the metric (HD 
represents Minkowski universe since the only non-trivial components of the curvature tensor are 
proportional to f^. Thus the simplest case for which (||) describes gravitational waves arises for 
/ = h(u)( 2 , where an arbitrary function h(u) characterizes the 'profile' of the wave. Solutions 
of this type are called plane waves (or homogeneous pp -waves) and have been investigated 
extensively in the literature (see 0,0 for Refs). This important class of exact radiative 
spacetimes has also been used for the construction of sandwich waves (see e.g. @H0) f° r 
which h(u) is non-vanishing on a finite interval of u only. Also, the very interesting problem 
of the collision of two plane waves has been thoroughly studied (see || for review of the topic 
and Refs.) 

Surprisingly, more general non-homogeneous pp -waves with generic or sandwich wave- 
profiles have not been investigated as much (a physical interpretation of some solutions of 
this type was proposed e.g. in [P],p!0|). On the other hand, much attention has been payed 
to impulsive waves described by the metric ([!]) with the function / of the form / = 5(u)F((), 
where S(u) is the Dirac distribution. These solutions have attracted attention recently not 
only for the investigation of gravitational radiation emitted during the collision of black holes 
|TTJ,||12[, but also in the context of field and string theories in a study of the scattering processes 



at extremely high Planckian energies, see [J3| and elsewhere. Such impulsive waves can natu- 
rally be understood as distributional limits of suitable sequences of sandwich waves regularizing 
the Dirac 5. Spacetimes of this type can also be obtained by boosting the Schwarzschild, other 
Kerr- Newman, and dilaton black holes [P^]- ||19|| , or axially symmetric solutions of the Weyl 
class [EDI to the speed of light. Some of these solutions can be interpreted as impulsive pp- 



waves generated by null particles of an arbitrary multipole structure . Another method for 
constructing general impulsive pp-w&ve spacetimes has been proposed in the 'scissors- 



and-paste' approach is based on the removal of a null hyperplane from Minkowski spacetime 
and re-attaching the parts by making a formal identification with a 'warp'. This corresponds 
to a specific coordinate shift |23| . 



Impulsive pp -waves can alternatively be introduced using a different coordinate system 
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in which the metric is explicitly continuous |TTJ, [P2"]-P5|. It has only recently been shown 
pj] that these approaches are equivalent in a rigorous sense. They lead to identical (unique) 
particle motion, i.e. the corresponding geodesies in ([!]), see ]2"7|]-[P5|], agree with those obtained 
from the continuous form of the metric. 

In our recent works |TT| we demonstrated that geodesic motion in the class of non- 



homogeneous vacuum pp-w&ve spacetimes is chaotic. This is the first example of chaos in 
exact radiative spacetimes (chaotic behavior of geodesies in various black-hole spacetimes is 
well known, see e.g. [p^l-[p6|, and references therein). It is the purpose of the presented 
paper to generalize these results which remained restricted to spacetimes (H) with / = - C( n , 
n = 3, 4, ■ • •, i.e. those having a constant profile, h(u) = C = const. In particular, we wish 
to investigate the behavior of geodesies in non- homogeneous sandwich pp -waves and in the 
corresponding impulsive limits. 

In section 2 we briefly summarize previous results concerning chaotic motion in pp -waves 
with a constant profile. Also, we introduce the concepts and quantities necessary for our 
subsequent investigation. Section 3 is devoted to geodesies in the simplest sandwich pp -waves 
having a 'square' profile. A rigorous characterization of the level of the fractal structure of basin 
boundaries (which separate initial conditions leading to different outcomes) by the 'number of 
radial bounces' of geodesies is introduced. This enables us to describe quantitatively how chaos 
smears as the sequence of sandwich waves approaches the impulsive limit. The same behavior 
of geodesies is investigated within the class of smooth sandwich waves in section 4, where a 
relation between the smearing of chaotic behavior and the gradual vanishing of the outcome 
channels is also discussed. In section 5 explicit (i.e. non-chaotic) geodesies in the resulting non- 
homogeneous impulsive pp- waves are presented and interpreted physically. Finally, in section 6 
we discuss the effect of the focusing of geodesies and describe the deformation of a disc of test 
particles. 

2 Chaotic motion in non- homogeneous pp -waves with a 
constant profile 



As shown in P0 |, |H] the geodesic equations for (|l|) reduce to 



C + ^C^ ' (2) 
u{t) = Ut + U , (3) 

v(r) = \U- X J [ 2CC - (/ + I) U 2 - e } dr + V , (4) 

where r is an affine parameter, dot denotes d/dr, U,U,V are constants, and e = —1,0,-1-1 
for timelike, null or spacelike geodesies, respectively. It suffices to find solutions of Eq. (0). 
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Introducing real coordinates x and y by ( = x + iy, this system follows from the Hamiltonian 

H= l(pl + pl)+V{x,y,u) , (5) 

where the potential is V(x, y,u) = \U 2 TZef. For non-homogeneous pp-waves given by 

f=lh(u)CT, n = 3,4,..., (6) 

the corresponding polynomial potential 

^^l^f/^W^eC (7) 

at any w = m is called 'n-saddle'. It can by visualized in polar coordinates p, <fi where 
( = pexp(i0), in which it takes the form V(p, </>, uq) = \U 2 h{uo) p n cos(n0). 

It was shown in a series of mathematical papers |3^|-P0|| that motion in the Hamiltonian 
system (|) with a polynomial potential (|7[) where h(u) = C = const, is chaotic. Note that, 
interestingly, in the simplest case n = 3 the corresponding 'monkey saddle' potential (after 
removing the factor CU 2 ) is V(x, y) = ^x 3 —xy 2 , so that we get the particular case of the famous 
Henon-Heiles Hamiltonian [[41|] (with missing quadratic terms) which is a 'textbook' example 
of a chaotic system. This was investigated by Rod |37j who concentrated on the topology of 
bounded orbits in the energy manifolds H = E > 0. The sets of orbits asymptotic to the basic 
unstable periodic orbits (denoted by H,) as r — > ±oo intersect transversely. This proves the 
existence of homoclinic and heteroclinic orbits and indicates the complicated structure of the 
flow. These results were later refined in [^] by showing that H,- are hyperbolic, so that they 
admit stable and unstable asymptotic manifolds. Finally, in |Mj the above Hamiltonian was 



presented as an example of a system for which the Smale horseshoe map can explicitly be 
embedded as a subsystem along the homoclinic and heteroclinic orbits. It was shown in |37] - 



40 that similar results hold for a general potential (0) with h = C . Therefore, geodesic motion 



in all non- homogeneous pp-w&ve spacetimes with the corresponding function (|6|) is chaotic. 

In order to support these arguments for the chaotic behavior, we investigated [[3(J-[[n|] the 
structure of motion by a fractal method. Complementary to the analysis described above, we 



concentrated on unbounded geodesies. The fractal method (advanced in relativity in |33j|- [|36[l 
and elsewhere) starts with a definition of several distinct asymptotic outcomes (given here by 
'types of ends' of all trajectories). Subsequently, a set of initial conditions is evolved until one 
of the outcome states is reached. Chaos is established if the basin boundaries which separate 
initial conditions leading to different outcomes are fractal. Such fractal partitions are the 
result of chaotic dynamics and measure an extremally sensitive dependence of the evolution 
on the choice of initial conditions. We have observed exactly these structures in the system 
studied. We integrated numerically the equations of motion given by (@), (f| for h = C . The 
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initial conditions (without loss of generality) were chosen such that the geodesies started (from 
rest) at r = from a unit circle in the (x, y)-plane. The initial positions were parametrized 
by an angle G [— 7r,7r) such that x(0) = cos0, y(0) = sin (p. In Fig. 1 we present typical 
trajectories for n = 3,4,5. Each unbounded geodesic escapes to infinity where the curvature 
singularity is located only along one of the n distinct outcome channels in the potential with 
the radial axis <ftj = (2j — l)7r/n, j = l,---,n (in fact, it oscillates around the axis with 
frequency growing to infinity and amplitude approaching zero |]30[]). These channels represent 
possible outcomes of our system and we label them by the corresponding values of j. In certain 
regions the function representing a portrait of the basin structure depends sensitively on 
initial position given by <fi — see Fig. 1 where the (finite resolution) results for n = 3, 4, 5 are 
shown. In the same diagrams we plot the function t s {4>) which takes the value of r when the 
singularity is reached by a given geodesic. The boundaries between the outcomes are fractal 
which we confirmed in on the enlarged detail, on the detail of the detail etc. up to the 
sixth level (where numerical errors became significant). At each level the structure has the 
property that between two connected sets of geodesies with channels ji and ji ^ ji there is 
always a smaller set of geodesies with channel j% such that j$ 7^ j\ and j'3 7^ ji- This has a 
counterpart in r s (</>), see Fig. 1. The value of r s diverges on each discontinuity of j (</>), i.e., 
on any fractal basin boundary. There is an infinite number of peaks, each corresponding to 
an unstable trapped orbit which never 'decides' on a particular outcome to infinity. Also, r s 
increases as one zooms into the higher levels of the fractal. This is natural since higher levels 
are generated by geodesies which undergo 'more bounces' in the inner region before escaping 
through one of the outcome channels. 



3 Motion in shock waves and smearing of chaos in the 
impulsive limit of simplest sandwich waves 

It can immediately be observed that the above results can easily be applied to a description 
of geodesies in shock pp -waves given by the metric (|]), (^) with h(u) = C0(w), where 0(m) 
is the Heaviside step function. It is natural to consider free test particles which are at rest 
(x — = y) in the flat Minkowski half-space u < 0. (The Minkowski coordinates are given by 
xm = V%x, Vm = V%y, zm = {v — u)/ y/2, tM = (v + u)/ v^2-) At u = these particles are hit 
by the shock and subsequently for u > they move in the wave with a constant profile h = C. 
The behavior of the corresponding geodesies has been summarized in section 2. In particular, 
the trajectories and the fractal structure of basin boundaries have again the form indicated 
in Fig. 1. Therefore, the geodesic motion in non-homogeneous shock pp-waves is chaotic in a 
rigorous sense. 

However, it could be argued that these above results concern very specific and rather 
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'unrealistic' classes of pp -wave solutions for which the profile function h(u) in (||) is constant 



on an infinite interval of the retarded time u. Such waves have an 'infinite duration' and 
a constant 'strength'. It is the purpose of the presented work to investigate more realistic 
non-homogeneous pp -waves, namely sandwich waves of this type described by functions h(u) 
having only finite support. 

For an investigation of sandwich waves it is convenient to parametrize the geodesies by the 
coordinate u instead of the parameter r; for ((u) we get from Eqs. @, @, (H) 



where prime denotes d/du. For timelike geodesies we can simply substitute u = Ur + U in the 
result in order to obtain the dependence on the proper time. 

It is natural to start with the simplest sandwich waves having a 'square' profile 



where a and a are positive constants. There are flat Minkowski regions in front of the wave 
(u < 0) and behind the wave (u > a). Within the wavezone (0 < u < a) the amplitude is 
constant, h = 1/5, and we can use the results described above. In particular, we can study 
a deformation of a ring of particles in the (x, ?/)-plane which are at rest in front of the wave. 
The particles start moving at u = and then follow exactly the same trajectories ((u) as 
shown in Fig. 1 for u < a. At u = a the influence of the sandwich wave ends, the potential 
(|7|) defining the outcome channels vanishes, and for u > a, the particles move uniformly in 
different directions along straight lines in the flat space behind the wave. Consequently, the 
fractal structure of basin boundaries indicating chaos is 'cut' at some level given by the value 
of a: there is not enough time for particles to bounce arbitrarily many times between the 
potential walls. As a —>■ 0, the number of bounces tends to zero so that the fractal structure is 
completely erased. In other words, for smaller a, the dependence on initial conditions is less 
sensitive since the prediction of the outcome can be done with only a lower resolution. The 
narrower the sandwich wave, the 'less chaotic' the corresponding geodesic motion. Finally, 
in the impulsive case given by the limit a — a — > of (^|), i.e. / = ^5(u)( n , n > 3, the 
motion is non-chaotic. This effect can be called a 'smearing of chaos' in the impulsive limit of 
non- homogeneous sandwich pp -waves. 

This behavior can be described formally. Let us define the number of bounces N of the 



geodesic as the number of (local) maxima of the function p(u) = \((u)\ = \Jx 2 {u) +y 2 {u) 
measuring the radial distance from the origin of the (x, |/)-plane. iV represents the number 
of times the geodesic crosses the phase-space surface of section p' = (with p" < 0) before 
escaping to infinity. This is a more appropriate measure here than the twist number used in 



C" + Ku)C 



0, 



(8) 




(9) 
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literature on chaotic dynamics []. As in the previous section, we consider a family of geodesies 
starting from rest from a unit circle p(0) = 1 with the angle <fi parametrizing their initial 
positions. In Fig. 2 we plot the function N(<p) for the geodesies given by Eq. (^j) with h = 1 
and n = 3. The sequence of graphs shows the zooming in of the fractal interval around the 
value ~ 0. Also, we plot the function u s {4>) which takes the value of u when the singularity 
at p = oo is reached by a given geodesic. It is obvious that the fractal structure described by 
the function iV(0) corresponds to the structure given by u s ((p). Also, N is a precise definition 
of the 'level' of the fractal since we have demonstrated that the iV-th level of the fractal 
structure is given by those geodesies which bounce exactly N times before they choose one of 
the outcome channels. Now we can describe how the fractal structure indicating chaos arises. 
Let us denote by Ik the set of initial conditions <fi generating geodesies with N(<f>) < k, where 
k = 0, 1,2, Clearly, J C I\ C Ii C • • • C (0,27r), and the complement of lim^oo Ik 
describes the fractal boundary intimately related to the fractal structure of the outcome basin 
boundaries (although the discontinuities in N((f>) and j(0) do not coincide). Now, we define a 
sequence of real numbers = max{u s (<f)) \ <fi G Ik}- Their physical meaning is the following: 
at u > at all geodesies generated by Ik have already fallen into the singularity at p = oo 
along the three channels given by <frj (after performing at most k bounces). From numerical 
simulations visualized in Fig. 2 we obtained the values do = 6.3, a\ = 9.3, <22 = 12.2, CI3 = 15.1, 
a 4 = 18.0, a 5 = 21.0. 

Now, we quantitatively characterize the behavior of geodesies in the sandwich wave given 

by (H). The amplitude within the wavezone is h = 1/a but it follows from Eq. (||) that we 

can obtain the corresponding values of u s ((p) from those for h = 1 (discussed above) by a 

simple rescaling u s (<p) = a/5m^ =1 (0), where w^ =1 (0) is drawn in Fig. 2. Notice that u s — > as 

h — y 00. It is obvious that all geodesies for which u s ((f)) < a, i.e. w(? =1 (0) < a/y/d,, have time 

enough to fall into the singularity (in the impulsive limit, a = a — > 0, this condition cannot 

be satisfied so that none of the geodesies considered ends in the singularity). Since iV(</>) is 

independent of the rescaling of u, the fractal structure of initial conditions given by Ik and the 

values of do not depend on a. Consequently, the condition a/yjli > ak guarantees that all 

geodesies starting from the initial set Ik will manage to fall into the singularity at p = 00 along 

the three channels (bouncing at most /c-times before that). For a sequence of sandwich waves 

with a = a leading to the impulsive limit as a — > 0, this condition reduces to a simple relation 

a > a\ which can be interpreted as follows. In order to emerge the fc-th level of the fractal 

structure of motion the sandwich wave must have at least the duration a such that a > a\. 

For example, if a > 441 then all geodesies with N < 5 bounces have time enough to choose the 

^^Note that actually the particle does not stop at this section — only its radial velocity vanishes — so that 
the event could better be called a 'turn' or a 'radial deflection' rather then a 'bounce' on the potential wall, 
but we use this last word for its natural intuitive meaning. 
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corresponding channel and fall along this into the singularity, i.e. J5 is fully developed. For 
smaller a the fractal structure is 'cut' at lower levels Ik so that the geodesic motion becomes 
'less chaotic'. For a < Oq ~ 40 the fractal structure vanishes since there is no time for any 
geodesic (performing at least one bounce) to choose the outcome channel and fall into the 
singularity. Thus, for a — > the motion is regular. This effect occurs in all non-homogeneous 
sandwich pp-waves with arbitrary exponent n = 3,4, • • -, only the specific values of for a 
given n are different: for higher n the values of are smaller. 

4 Motion in smooth asymptotic sandwich waves 

We have observed that, in non- homogeneous sandwich pp -waves with a short duration, many 
particles still remain in the inner region when the potential (]?]) defining the outcome channels 
vanishes. These particles move through various points in different directions having different 
velocities when the sandwich wave ends. Subsequently, they move uniformly in the Minkowski 
space behind the wave and it is obvious that they will not follow former channels to the 
singularity at p — 00. In fact, any point is accessible by some trajectory. The narrower the 
sandwich wave, the greater the number of particles moving 'outside' the channels. 

In order to illustrate such behavior of geodesies (and to emphasize this important aspect of 
the impulsive limit) we consider another type of sandwich waves (||) with the profile function 



where b, b are positive constants. Note that the corresponding radiative spacetimes are curved 
everywhere, becoming flat (exponentially fast) only asymptotically as u —>■ ±00. On the other 
hand, these waves are smooth since all derivatives of h(u) are continuous. Moreover, (|i~0"D is 
more appropriate for numerical integration of geodesies (contrary to the case of simplest waves 
given by Eq. @ we need not apply the junction conditions which guarantee the continuity of 
motion at u = and u = a where h(u) is discontinuous). Again, we consider a ring of free 
test particles which are at rest at u = (other initial conditions, e.g. a ring which is at rest at 
u = —00 lead to analogous results). In Fig. 3 we show their geodesic motion in the wave with 
n = 3 and h(u) of the form (0), for different values of the parameter b = b. For very small 
values of b the profile remains almost constant on a large interval of u and the trajectories 
'coincide' with those shown in Fig. 1. The outcome channels are very well defined and the 
motion is chaotic. However, with a growing value of b these channels become fuzzy. It is 
harder and harder to distinguish between them, so that the basin boundaries lose their fractal 
structure. For b — > 00, the profile function flip]) approaches the Dirac delta (in a distributional 
sense) and again, in the impulsive limit the geodesic motion is not chaotic but regular. 




(10) 
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This effect is quantitatively characterized in Fig. 4 where we plot the function j(0) for 
different values of b. For this purpose we define three outcome windows as small intervals A0 
of angles around the radial axes 4>j of the three outcome channels (localized sufficiently far away 
from the origin); here we consider \<f> — <f>j\ < A0 = 0.1. When a geodesic starting at <fi passes 
through some outcome window we assign the corresponding value of G (1, 2, 3) to it. If the 
trajectory does not pass through any of the three windows then we define j (<j)) = which means 
that the geodesic does not approach infinity along the outcome channels <pj. From a sequence 
presented in Fig. 4 it is obvious that the number of geodesies localized outside the channels 
increases with a growing value of b, i.e. as the sandwich profiles approach the impulsive limit. 
In the same diagrams we also plot the function u s (<p). One can observe that the number of 
geodesies reaching the singularity in finite values of u s rapidly decreases with a growing b. In 
particular, all geodesies outside the channel windows (with j(<f)) = 0) have u s (<f)) = oo. There 
is a critical value b c ~ 0.09 such that for b > b c all geodesies considered have u s = oo, see 
Fig. 4. These geodesies reaching p = oo at u s = oo are, in fact, given by specific solutions of 
Eq. (|8]) for which the condition |ft.(u)l <C p x ~ n is satisfied for large values of u. Consequently, 
< |C"I ^ 1) \u"\ ^ IC"I ^ I; i- e - the acceleration is negligible and asymptotically the 
particles move uniformly along straight lines. The coordinates x and y depend linearly on u, 
so that p = oo is reached at u = oo. Note that for the smooth profile ( |T0"D of the wave with 



n = 3 this condition reduces to cosh(26n) ^> \Jb p{u). Once this condition is satisfied for some 
u, it is valid for all greater values of u since the left hand side grows exponentially while the 
right hand side only linearly as u — > oo. 

Finally, we wish to find the critical value of b c . From Fig. 4 it is obvious that for b = b c 
even the geodesic starting at u = from rest from the unit circle at = tt (and similarly from 
= ±7r/3) will reach p = oo with an infinite value of u s . In real coordinates, this geodesic is 
given by x" + h{u)x 2 = 0, y(u) = 0, with x(0) = —1 and x'(0) = 0. Introducing ip — In \x\ 
we get the equation ip" = h(u)exp(ip) — ip' 2 , ip(0) = = ip'{0); the corresponding solutions 
for various values of b = b are plotted in Fig. 5. For small b the function ip(u) is convex 
and diverges at finite values of u. For large values of the parameter b, it becomes concave 
so that the geodesies approaches infinity only as u — > oo. The boundary between these two 
types of behavior defines b c . From numerical simulations we obtained an approximate value 
b c = 0.0872374, see Fig. 5. 

5 Motion in impulsive waves 

We have demonstrated above that chaos disappears when the sequence of non-homogeneous 
sandwich pp -waves approaches the impulsive limit. In this section we concentrate on a descrip- 
tion of motion in these impulsive spacetimes. In particular, we present and discuss an analytic 
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solution to the geodesic equations, thus proving explicitly that the motion is fully non-chaotic. 

For this purpose it is convenient to use a coordinate system for impulsive pp-w&ve space- 
times which is continuous for all values of u [2~4"|, [25|. Following [ 25fl , the metric can be written 



in the form 

ds 2 = 2 |df7 - \uQ{u) [d 2 F{ V )/dr] 2 ] dr/| 2 — 2 dudr . (11) 
The transformation relating (|ll]) and (|l|) with / = 5(u)F(() is 

C = V~ \uQ{u) (dF/dfj) , (12) 
v = r - + F) + \uQ{u) \dF/dri\ 2 . (13) 

Obviously, there are privileged geodesies rj — rjo =const. in ([11]) corresponding to free particles 
which remain at rest in flat Minkowski half-space u < in front of the impulsive wave. After 
the passage of the impulse these geodesies are still given by rj = r/o but the flat half-space 
u > behind the wave is naturally described by coordinates £ and v (see ([!])) in which the 
motion is uniform. From Eq. ( |12"D it follows that, for impulsive waves given by F(() = ^C( n 
(so that h(u) = C5(u)), the motion in the transversal plane is described simply by ((u) = 
r/o — CQ{u) uviq' 1 . If we parametrize the initial position of each particle by ((u < 0) = t]q = 
p o exp(i0 o ), the motion behind the impulse {u > 0) is explicitly given by 

x{u) = Pq cos 4> — Cpg" 1 cos[(n — 1)0q] u , 

y(u) = p sin O + CX™ 1 sin[(n - 1)0 O ] u , (14) 

so that 

p 2 (u) = p 2 + C 2 pl n - 2 u 2 - ICpl cos(n0 o ) u . (15) 

The motion is uniform, i.e. the trajectories are straight lines with the velocity of each particle 
(x',y') = Cpo _1 (— cos[(n — l)0o]>sin[(n — l)0o]) being constant. Thus, the inclination of the 
straight trajectory in the transversal (x, y)-plane is a = tan(y' / x') = tt + (1 — n)<p , and the 
speed is a/ x' 2 + y' 2 = \C\ Pq^ 1 . Notice that the speed depends on p while the direction of 
motion on <f) only. The above geodesies can easily be visualized. We do not present their 
trajectories which would be very similar to those in sufficiently narrow sandwich waves, such 
as the b = b = 0.5, n = 3 case shown in Fig. 3. Instead, we draw in Fig. 6 'complementary' 
pictures showing a deformation of a ring of free test particles (for po = 1 = C, n = 3). The 
initial circle at u = is continuously deformed into smooth curves with n growing loops which 
are visualized here as sequences of 4, 8, 13 and 22 consecutive steps Am = 0.1, i.e. the largest 
connected curve describes the deformation of the ring at u — 2.2. Similarly, in Fig. 7 the 
deformation is shown for n = 3,4,5,6 and u = 0.25,0.5,0.75,1. The most distant particles 
are those which started at 0o = 4>j — (2j — l)n/n, j — 1, 2, • • • , n with cos(n0o) = —1 so that 
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p{u) = \po + CpQ~ 1 u\; they exactly follow the radial axes of the outcome channels which would 
be present for corresponding pp-w&ves with a constant profile. On the other hand, for 0o given 
by cos(n</>o) = 1, we get p(u) = \po — Cp^~ u\. 

It can be observed that for small values of u the deformation of the ring agrees with that 
shown in Fig 8. of Ref. |TT| describing chaotic motion in the non- homogeneous pp -waves with 



a constant profile. The principal difference occurs for large u. In the impulsive case there are 
no subsequent loops arising from more and more bounces in the inner region so that the circle 
can not be deformed in a fractal way with different segments moving to different outcome 
channels. Instead, there are no channels to the singularity, the trajectory of each particle is 
explicitly given by ([14]), and the motion is non-chaotic. 

These results can easily be generalized to the case when the test particles are not at rest 
initially. From the metric ([!]) it is obvious that, using the coordinate (, the motion must 
be uniform in both flat half-spaces, i.e. ((u < 0) = ip u + r] , ((u > 0) = Xo u + Co> where 
■00; Voi Xo, Co are constants. The relation between these parameters can again be found using 
the continuous form of the impulsive metric (0). By solving the corresponding geodesic 
equations we could obtain f](u). Although this function is very complicated for u > 0, it 
has the property that rj(u) and rj'(u) are continuous (even at u — 0). Using this fact and 
Eq. (0), which for F(() = -C( n reduces to ((u > 0) = rj{u) — Cur] n 1 (u), we obtain 



lim u ^ C(w > 0) = Co = 77(0) = r/ and {'(u > 0) = X o = rf(0) - Crf^W = ^o - Cfj^~\ i.e. 

au) = VPo-CQ{u)f 1 r 1 }u + r lo . (16) 

Therefore, the trajectory of each particle in the transversal plane is a continuous but refracted 
straight line with a discontinuity in the velocity at u = given by A(' = xo ~ ipo — ~C fjQ -1 . 
Notice that the value of the jump depends on the position rj of the particle at u = only, 
not on its actual velocity. By parametrizing rjo = po ex P(i0o) — Co we immediately obtain 
the discontinuity in velocity, (Ax', Ay') = Cp r ^ l {— cos[(n — l)</> ], sin[(n — 1)0q]). Thus, the 
result is trivial: if the particle is not at rest initially, behind the wave its constant velocity 
merely superimposes to the effect of a characteristic jump in velocity given by the impulse. 
By denoting cota^ = x'(u < 0) = TZeipo, cot (3 X = —x'(u > 0) = — TZexo, an d similarly 
cota^ = y'(u < 0) = Imipo, cot (3 y = —y'{u > 0) = —Imxo, we can rewrite the expression 
for AC' as cot a x + cot (3 X = Cp$~ cos[(n — l)0o] an d cot a y + cot (3 y = —Cp$~ l sin[(n — l)0o]), 
which generalize to non-axisymmetric cases the 'refraction formula' for deflection on (null) 
geodesies in the axisymmetric Aichelburg & Sexl spacetime [0 (see e.g. P5[ , |42| 



So far, we have concentrated on a description of motion in the transversal plane. We should 



also comment on behavior in the longitudinal direction. From the form of the metric (|TT|) it 



follows that, for geodesies T] = r] , the coordinate r(u) is given by r{u) = s u + r , s ,r being 



11 



constants. Moreover, r(u) is continuous for all u. Due to the Eq. (0), there is a jump in 



v(u) = s u + r + C 2 p 2 n ' 2 uQ(u) - ^Cp^ cos(n0 o ) &(u) (17) 

at it = given by Av = — yCp 1 ^ cos(n0o) depending both on p and 4>o- Since the Minkowski 
coordinates are zm(u) = -^[v(u) — u], t M {u) = ^[v(u) + it], there is a discontinuity Azm = 
AtM = ~^^v at u = 0. This 'shift' effect of impulsive wave on geodesies is well-known for the 
axisymmetric Aichelburg & Sexl spacetime flliU , see [fUH , p3| , ||42| , and elsewhere. 



A solution of the geodesic (and geodesic deviation) equations in general impulsive pp- 
wave spacetime has been presented in p7]-||29||. Starting from the distributional form of the 



metric (|l|), one has to deal with ill-defined products of distributions. In order to obtain 
correct results in a mathematically rigorous fashion, careful regularization procedures and 
delicate manipulations during calculation of the distributional limit are required. In fact, the 
Colombeau theory of generalized functions providing a suitable consistent framework has to 



be applied. Using this rigorous solution concept, it was shown in |29| that in the impulsive 
limit the geodesies are totally independent of the regularization, i.e. on the particular shape 
of the sandwich wave (the impulsive wave 'totally forgets its seed'). Here we obtained an 
identical explicit form of geodesies starting from the continuous form of the impulsive metric 



(|TT|). It was demonstrated in |26[ that these two approaches are equivalent in a mathematical 
sense (even if the transformation (112]), (pf) is discontinuous). However, our main goal here 
was not to re-derive the known geodesies but to concentrate on their physical description and 
visualization. 

6 Focusing of geodesies and caustic properties 

Obviously, geodesies with parallel trajectories in flat half-space in front of a sandwich or impulse 
wave are refracted. A natural question arises whether a specific character of the refraction leads 
to some form of focusing of the corresponding geodesies. Such an astigmatic focusing effect is 
well-known for plane waves, and results in interesting caustic properties thoroughly investigated 
in 0. For principal reasons (non-existence of a diagonal Rosen form of the metric, see Eq. 
(15) in p6| , or the presence of chaos in geodesic motion) it would be a very difficult task to 
reproduce these results for non-homogeneous sandwich waves. Therefore, we restrict ourselves 
to impulsive gravitational waves only. 

It immediately follows from Eqs. (|H|), (|I7|) that particles staying at fixed xo, yo, z in front 
of the impulse, move at u > according to 

x M (u) = V2x(u) , Vm(u) = V2y(u) , (18) 
V2 [z M (u) - Zq] = C 2 p 2 n - 2 u - \Cpl cos(n0 o ) , (19) 
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where xm,Um, z m are Minkowski coordinates behind the impulse, x(u) and y{u) are given 
by (|14|). In particular, for impulsive plane waves {n = 2) we get xm(u) = (1 — Cu)xq, 
Vm(u) = (1+Cu)y and 2\/2 [z M (u) — z ] = C(l+Cu)y 2 — C(l— Cu)xl, where x = y^PoOos 4> , 
yo = v^posin^o- Consequently, at u = Uf = 1/C > one gets xm(w/) = 0, ?/m(m/) = 2y 
and v2[zjvf(^/) — ^o] = CVo> so that all particles meet at xm = 0, independently of their 
initial position xq in the ^o-plane. The generically quadratic surface (hyperbolic paraboloid) 
u =const. degenerate at u = Uf to a parabolic caustic line % = 0, % = jj0m + ^o- This 
is shown in Fig. 8 which visualizes the coordinate singularity formation. It plays a crucial role 
in colliding plane-wave solutions where the corresponding global 'fold' singularities arise, as 
described in Consequently, these spacetimes are globally hyperbolic in contrast to single 
plane-wave solutions where the focusing effect on null cones forbids an existence of a spacelike 
hypersurface which would be adequate for the global specification of Cauchy data |15[ . 



Here we wish, however, to investigate key features of the focusing effect in more general 
non-homogeneous impulsive pp-wave spacetimes with n = 3,4,5, •• •. From the deformation 
of a ring of free particles indicated in Fig. 7 it follows that only n privileged particles can 
collide at the origin (in contrast to the plane-wave case where particles with arbitrary xq are 
focused along the xm = line). It can easily be derived from Eq. ( |1"5|) by setting p{u) = 
that only those n particles with cos(n0 o ) = 1 reach the origin of the transversal (x, y)-plane, 
simultaneously at u — Uf = p^~ n /C. Notice that particles starting from a larger circle focus 
sooner. This is understandable: they have to travel a bigger distance po, however, their speed 
is also bigger, Cpg -1 , so that the travel time is indeed Po~ n jG. Note that the occurence of such 
focusing is associated with the coordinate singularity in the Rosen form of the metric (|TT|). 

At u = u c the circle is deformed into a curve with n cusps, see Fig. 7. The value of u c 
can be found from the condition (<90/<90 o )|cos(n0o)=-i = 0> where = aicta.n[y(u) /x(u)}, with 
x(u),y(u) being given by (^4|). This leads to equation (n—l)C 2 pl n ~ 2 u 2 c + (n — 2)C ' p^u c — p\ = 
admitting a unique positive solution u c = p1~ n /[C(n — 1)] = Uf/{n — 1). For p = 1 = C we 
get a simple expression u c = l/{n — 1) which is in agreement with Fig. 7. 

In fact, the shape of deformation of the initial circle p = po of test particles in the transversal 

plane has a beautiful geometrical interpretation: at any time u, the circle is deformed into a 

curve called a hypotrochoid. Such a curve is generated by a point p attached to a small circle of 

radius B rolling around the inside of a large fixed circle of radius A, with H being the distance 

from p to the centre of the rolling circle P6|. The Eqs. ( ]I4]) for particle motion behind the 

impulsive wave are just the parametric equations for a hypotrochoid with the identification 

n 1 

A = -po, B = -po, H = Cup n Q -\ (20) 

n — 1 n — 1 

where 0o is the parameter. We observe that A = nB so that the hypotrochoids are closed 

curves with n loops in our case. Also, the parameter H grows linearly with u. At u = we 
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have H = and the curve is just the initial circle of radius A — B = p . At u = u c the 
parameter has the value H = p /(n — 1) = B and the hypotrochoid reduces to the curve called 
a hypocycloid which has n cusps. Note that it degenerates for n = 2 to a line segment, a 
3-cusped hypocycloid is called a deltoid, a 4-cusped hypocycloid is called an astroid (see Fig. 7 
for n = 3 and n = 5 case). Finally, at the focusing time u = Uf > u c we get H = po = A — B 
so that the hypotrochoid n-times intersects the origin. The corresponding curve is called a 
rose since it resembles a flower with n petals (for n = 2 the rose degenerates to a line since 



It should be emphasized that the above behavior of geodesies in the (x, y)-plane describes 
only part of the overall deformation. We have to consider not only the transversal deformation 
but also the motion in the ^/-direction described by (]T9|) , as we have done for the plane- wave 
case in Fig. 8. At u = the longitudinal deformation \/2[zm(u) — Zq\ caused by the impulse 
is given by the second term on the right hand side of Eq. (fl9|), which is exactly the shift 
Av = — fCpo cos(n0 o )- For u > the first term increases with u; at u w = j l Po~ n /C = ^-Uf the 
deformation of the whole disc p < p , z =const., becomes non-negative and subsequently grows 
to positive values, uniformly with u. In combination with the deformation in the transversal 
plane we get an interesting new effect of 'wrapping up' of the inner part of the disc into a closed 
caustic surface as shown in Fig. 9. Note that this behavior is specific for non-homogeneous 
pp -waves and does not occur in impulsive plane waves since u w = Uf when n = 2. 

7 Conclusions 

We have shown for non-homogeneous pp-w&ve spacetimes that the motion of free test particles 
is chaotic. However, for non-homogeneous sandwich gravitational waves the chaotic behavior 
smears as the duration of the wave tends to zero. In the limit of impulsive waves the motion 
is regular. The focusing effect of non-homogeneous impulsive waves can also be described ex- 
plicitly and leads to a formation of a specific closed caustic surface. It would be an interesting 
task to investigate from this point of view geodesies in other classes of exact radiative space- 
times and elucidate an open question whether such a type of behavior is common in general 
relativity. 

Also, it is obvious that sandwich pp-w&ve solutions may serve as local models of gravi- 
tational waves far away from a radiating source. In realistic situations these waves should 
contain non-homogeneous components whose influence on test particles has been described 
above. In particular, chaotic-type effects on null geodesies could in principle be observable 
after the passage of the sandwich wave as characteristic "chaotic deformation patterns" of an 
observer's view of the sky, i.e. as peculiar changes of positions of stars and galaxies. These 
specific astronomical consequences will be investigated elsewhere. 
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Figure Captions: 



Fig.l Geodesies starting from a unit circle escape to infinity only along one of the n 
channels (left). The functions j(<f>) and T s (<f>) indicate that basin boundaries separating different 
outcomes are fractal (right). 

Fig. 2 Plots of the functions N(<f>) and u s (<f>) in the sequence of ever narrower intervals 
exhibit their fractal structure. Obviously the number N of bounces of each geodesic is an 
appropriate measure of the fractal level. 

Fig. 3 Geodesies in the smooth sandwich gravitational waves given by Eq. (10), n — 3. For 
small values of the parameter b = b, most of them escape to infinity only along the well-defined 
three channels. For larger values of b these channels become fuzzy and the basin boundaries 
separating different outcomes are losing their fractal structure. In the impulsive limit {b — > oo) 
the motion is regular. 

Fig. 4 Plots of the functions j((f>) and u s (<f>) describing behavior of geodesies in smooth 
sandwich gravitational waves. When the parameter b = b is small, most geodesies reach the 
singularity in finite value of u s . For b > b c m 0.0872374 all geodesies have u s = oo. Also, for 
higher b the number of geodesies outside the three outcome windows (j = 0) is greater. 

Fig. 5 The function ip(u) = ki\x(u)\ describing specific geodesies changes its character at 
the critical value b c xs 0.0872374. For b > b c the geodesies approaches p = oo only as u — > oo. 

Fig. 6 The deformation of a ring of free test particles which are at rest initially in the 
tranvsveral (x, y)-plane by the influence of an impulsive gravitational wave with / = |<5(w)C 3 at 
u = 0.4, 0.8, 1.3 and 2.2 with steps Au = 0.1. Bellow is the 3D visualization of the deformation 

(x(u),y(u)). 

Fig. 7 Deformation of a unit ring of particles by impulsive waves with / = ^5(u)( n at 
u = 0.25,0.5,0.75, 1 for n — 3,4,5,6. At u — Uf = pl~ n /C = 1 all the privileged n particles 
focus at the origin of the transversal plane. The curves are hypotrochoids, as explained further 
in the text. 

Fig. 8 Deformation of a disc consisting of free particles p < 1, z — 0. By an impulsive 
plane wave (n — 2, C — 1) the dies is deformed into a hyperbolic paraboloidal surface at 
u =const.> in the Minkowski space behind the wave. At u — Uf — 1 the surface degenerates 
to a caustic parabolic line. 

Fig. 9 Deformation of a disc p < 1, z = of particles by an impulsive no n- homogeneous 
wave / = |£(-u)C 3 at u — 0,0.25,0.5,0.75,1 and 1.25 in the Minkowski space behind the 
wave. At u — Uf — 1 all the privileged three particles collide and with the deformation in the 
longitudinal ^M-direction a closed caustic surface is created. 
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